Experimental and materials considerations for the topological superconducting state 
in electron and hole doped semiconductors: searching for non-Abelian Majorana 
modes in ID nanowires and 2D heterostructures 
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In proximity to an s-wave superconductor, a one- or two-dimensional, electron- or hole-doped 
semiconductor with a sizable spin-orbit coupling and a Zeeman splitting can support a topological 
superconducting (TS) state. The semiconductor TS state has Majorana fermions as localized zero- 
energy excitations at order parameter defects such as vortices and sample edges. Here we examine 
the effects of quenched disorder from the semiconductor surface on the stability of the TS state 
in both electron- and hole-doped semiconductors. By considering the interplay of broken time 
reversal symmetry (due to Zeeman splitting) and disorder we derive an expression for the disorder 
suppression of the superconducting quasiparticle gap in the TS state. We conclude that the effects of 
disorder can be minimized by increasing the ratio of the spin-orbit energy with the Zeeman splitting. 
By giving explicit numbers we show that a stable TS state is possible in both electron- and hole- 
doped semiconductors for experimentally realistic values of parameters. We discuss possible suitable 
semiconductor materials which should be the leading candidates for the Majorana search in solid 
state systems. 

PACS numbers: 03.67.Lx, 71.10.Pm, 74.45.-l-c 



I. INTRODUCTION 

In a recent paper— it has been shown that a 2D spin- 
orbit coupled electron-doped semiconductor, in proxim- 
ity to a bulk s-wave superconductor and an externally 
induced Zeeman splitting, can support a topological su- 
perconducting phase with Majorana fermion modes at 
vortex cores and sample edges i.e. at order parameter de- 
fect locations. It has also been realized from dimensional 
reduction that the edge Majorana modes in a sample with 
rectangular geometry (of width W) turn into localized 
end Majorana modes in the corresponding ID geometry 
[W — > 0) and, further, in ID the excitation gap above the 
Majorana modes (mini-gap) should scale as the induced 
superconducting gap A (thus avoiding the problem of a 
tiny mini-gap, A^/ej? ^ 0.1 mK, as in 2D (p^^ -|-ipj,)-wave 
superconductors). The end Majorana modes in the ID 
geometry should be experimentally visible in local zero- 
bias tunneling experiments. These results, numerically 
confirmed soon thereafter,^ have appeared with analyti- 
cal and numerical details for both ID and 2D semicon- 
ductors in Ref. Q- It has also been pointed ovA^ that 
the required Zeeman splitting in ID can be induced by 
a magnetic field parallel to the adjacent superconduc- 
tor (thus avoiding the problem of orbital effects^), and 
the non-Abelian character of the end Majorana modes 
can be probed by Josephson experiments. Very recently 
it has been shown that the generic Luttinger Hamilto- 
nian applicable to the hole- doped nanowires also sup- 
ports end Majorana modes - in a manner similar to its 
electron-doped counterpart. The ID electron- or hole- 
doped wires in the TS state can be arranged in a quasi- 
ID network geometrj«^to test non-Abelian statistics^ and 
perform topological quantum computation (TQC)^'^*^ in 



the Bravyi-Kitaev (BK) schemei^i In principle, TQC in 
the BK scheme is also possible with the 2D semiconduc- 
tor TS states using Majorana fermion interferometrj*i^ 
analogous to that in the v = hjl fractional quantum 
Hall (FQH) stateSfi^ chiral p-wave superconductors^^ 
and TS states on the surface of topological insulators J^ii^ 
Efforts to realize a semiconductor TS state with Majo- 
rana fermions in proximity to s-wave superconductors 
are currently underway in many laboratories world-wide, 
concentrating on both ID semiconducting nanowires and 
2D semiconductor heterostructures in close proximity to 
a regular bulk s-wave superconductor (e.g. Al, Nb). 

In view of the ongoing experimental efforts on semi- 
conductor TS states, it is important to understand the 
effects of disorder from the semiconductor surface on the 
realizability of the TS state. (In recent worksii it has 
been established that the disorder residing in the ad- 
jacent bulk superconductor has negligible effect on the 
semiconductor TS state, and therefore in this paper we 
will ignore the effects of disorder from the bulk super- 
conductor.) This is an especially important questioni^Ti^S 
because the TS states in semiconductors explicitly break 
the time reversal (TR) symmetry and it is known that 
the superconducting quasiparticle gap in such systems, 
unlike that in TR-invariant superconductors is sup- 
pressed by disorder. Below we will measure the extent 
of TR breaking in the semiconductor TS state (in both 
ID and 2D) by the ratio r = Vz I {ctkF)-, where Vz is the 
Zeeman energy in the TS state and akp gives the typi- 
cal spin-orbit energy scale (with a as the Rashba spin- 
orbit coupling constant). We will show that the disor- 
der suppression of the superconducting quasiparticle gap 
from its clean value increases with increasing values of r. 
Nonetheless, by giving explicit numbers we show that a 
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reasonable TS state quasiparticle gap can be experimen- 
tally achievable in both electron- and hole-doped systems 
even in the presence of realistic disorder. For electron- 
doped wires we find that with a mobility of ^ 100, 000 
cm^/V-s a measurable robust TS state gap of 50 — 100 
mK is achievable. We also find that for hole-doped wires 
with a significantly larger spin-orbit energy Eso ~ 30 
meV,'^'^'^'^ a larger gap of 0.8 K is achievable for a mobil- 
ity of - 100 cmVV-s. 

In what follows, by we denote the quasiparticle gap 
(pair potential) in the (proximity-inducing) bulk super- 
conductor. For the semiconductor (2D or ID), a is the 
Rashba spin-orbit coupling constant, m* is the effective 
mass (of electrons or holes), Vz is the Zeeman splitting 
required for the TS state, Aq is the proximity-induced 
pair potential, A is the quasiparticle gap at the Fermi 
surface, kp is the Fermi momentum {h = e = ks = 1), 
and 1/Tsm is the effective intrinsic semiconductor disor- 
der scattering rate. In the semiconductor, the combina- 
tion of the Zeeman splitting and the spin-orbit coupling 
results in a Fermi-surface quasiparticle gap A, which dif- 
fers from the proximity-induced pair potential Aq, and 
thus we have three distinct superconducting gaps As, A 
and Aq to consider in this problem. 

In Sec. II we first give our starting Hamiltonian for the 
semiconductor which includes a proximity-induced pair 
potential, an externally applied Zeeman splitting and a 
disorder potential. We then calculate the dependence of 
the disordered TS state quasiparticle gap on the semi- 
conductor disorder scattering rate 1/Tsm (Eq- HH])- In 
Sec. Ill, we review the two systems that we will apply 
our results to, namely, electron and hole-doped semicon- 
ductors. In section IV, we use the results of Sec. II and 
Sec. Ill to give explicit realistic values for the various 
parameters for a stable TS state in electron-doped semi- 
conductors. This is followed by the case of hole-doped 
semiconductors, where again using explicit values of the 
parameters we show that a stable TS state is possible. 
We find that the TS state gap in the hole-doped case can 
be much larger and more robust to disorder effects than 
in its electron-doped counterpart due to larger values of 
a and to*. We conclude in Sec. V with a summary and 
an outlook with a detailed qualitative discussion on why 
the key problem of choosing the optimal materials com- 
bination for realization of the TS state is experimentally 
important. 

II. BDG HAMILTONIAN WITH DISORDER 

The BdG Hamiltonian in real-space for spin-1/2 elec- 
trons with spin-singlet pairing potential is written as 

where Aq is a real, constant, s-wave spin-singlet pair- 
ing potential in the semiconductor which is proxim- 
ity induced from an adjacent superconductor. Ham = 



ffsm('", c; rV) is the non-interacting part of the semi- 
conductor Hamiltonian with disorder and tp^sm is the 
chemical potential. Here Ox,y,z are the Pauli matrices. 
For a disordered semiconductor, with spatially localised 
(i.e. short-ranged) spin-independent disorder, i/s,„ can 
be written as 

Hsm{r,a;r'a') = Hs,n,ciean{r,a]r'a')+V{r)5{r-r')Saa' 

(2) 

where V(r) is the local potential induced by disorder. 

In the limit of weak disorder and weak pairing where 
both V{r) and Aq are smaller than the inter-band spac- 
ing of the spin-orbit bands in the semiconductor, it is 
convenient to work in the basis of Bloch eigenstates |n, k) 
of Hsm,cieanifO';r'a'). In particular, in the case where 
a single band, say n, crosses the Fermi-surface (which 
is necessary for realizing a topological superconducting 
state2ii2^), we can calculate the low-energy spectrum of 
the Hamiltonian from an effective Hamiltonian which is 
obtained by taking matrix elements of the BdG Hamil- 
tonian in Eq. [1] with respect to the Nambu spinors 

l'^'fe>Wam6« = e\n,-k) ) 

where Q — iUyK with K being the complex conjuga- 
tion operator. In this approximation, the translationally 
invariant part of the BdG Hamiltonian has the form 

Ho(fc, fc') = ( "^Za*-'^ u f'-e. '^('^ - '^') 

(4) 

where we have suppressed the band-index n, since 
we are restricted to a single band and the term 
Hsm.,ciean{i'o-',f'o'') has bccn replaced by its eigenvalue 
Efc in the relevant band. In the electron-doped case, 
includes effects on the dispersion of the parameters Vz 
and a (see Eq. [H]) while in the hole-doped case, it addi- 
tionally contains information about the Luttinger param- 
eters (see Eq. [29|l. We note that, in spite of its apparent 
complexity involving the superconducting proximity ef- 
fect, spin-orbit coupling, Zeeman splitting, and disorder, 
we are still dealing here with an effective exactly solvable 
one-particle quantum problem. 

For weak superconducting pairing and disorder scat- 
tering, the TS quasiparticle gap is determined by the en- 
ergies and wave-functions near the Fermi energy (ip^sm- 
In particular, we will approximate the dispersion around 
the Fermi energy by ~ ^p.sm + WF(|fc| — kp) where vp 
is the Fermi-velocity. Moreover, since the relevant states 
are near the Fermi wave-vector (i.e |fc| « kp), we will 
assume that matrix elements such as (fe|fci) depend only 
on the directions of the momenta fc and fci . To calculate 
the matrix elements of |fc), we will assume that the Bloch 
eigenstates have a simple decomposition as 

(r; a\k) = Uk{r; a) = e''^ '' Uk{a) . (5) 

This assumption is valid as long as all perturbations in- 
volve momentum transfers that are smaller than a Bloch 
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vector. For electrons a is simply the electron-spin. When 
applied to holes a would represent the 4-component 
pseudo-spin degree of freedom associated with the Lut- 
tinger model.- 

In this paper we will consider both ID and 2D sys- 
tems, and these cases will turn out to be closely related. 
In the two-dimensional case, we will restrict ourselves to 
rotationally symmetric systems with circularly symmet- 
ric Fermi surfaces where k = [kp cos Ok^kpsuvOk)- Using 
the azimuthal symmetry of the Bloch Hamiltonian corre- 
sponding to Hg^ cieam we Can define a Hermitian matrix 
Rz, which generates z-axis rotations so that 

=e*(«^^^)(«'=)||fc|,0) (6) 

where A is chosen so that ||fc|, 0^) is a single valued func- 
tion of 6k. Since Kramer's theorem requires Rz to have 
half-integer eigenvalues A = 1/2. For one-dimensional 
systems, such as electron- or hole-doped semiconduct- 
ing nanowires, we will consider systems with inversion 
symmetric Fermi points ±kp. In this case, the relevant 
states at the Fermi-level are |±fci?). The 2D systems 
with Rashba spin-orbit coupling will turn out to have 
results that are solely dependent on the wave-functions 
only at \kF,9k = 0, vr), which are identical in form to the 
states IzLkp) relevant for the ID case. Therefore, these 
two cases are closely related and for results such as the 
quasiparticle gap in Eqs. [T7] and I18[ we will only state 
the result for the ID case and imply a similar result in 
2D, which can be derived easily. 

The superconducting pairing potential in the pro- 
jected Hamiltonian in Eq. 2] is the matrix element of Ao 
in Eq. [T] between states near the Fermi-level, which is 
written as 

Afc = Ao(fc|e|-fc). (7) 

Using Eq. [S] in the rotationally symmetric 2D case, Afc 
simplifies to 

Afc = Ao (fc| e \-k) = Ao (0| e-'(-R--^)«'«ee*(^--^)'''= [tt) 
= e^'^^-^Ao (0| e Itt) = Ao (0| 6 {tt) e*''^ (8) 

where 9^ is the angle of the wave- vector k in the k^ — ky 
plane. Here we have assumed that i?^, which transforms 
kp — > —kp, is odd under time-reversal symmetry. The 
TS quasiparticle gap for the Hamiltonian Eq. 21 in the 
absence of disorder (i.e. V{r) = 0), is given by 

A = |Afc|. (9) 

The disorder scattering, induced by the potential V{r) 
in Eq. [21 can be projected into states near the Fermi 
energy in a similar way as 

where v{k, k') = {k\ V{r) \k'). For spatially uncorrelated 
white noise disorder (i.e. {V{r)V{r')) = v'^5{r — r')), and 



using the simple form for the Bloch functions Eq. [5l the 
disorder propagator (u(fci, A;2)w(v3, ^4)) can be written 
as 

(V(fel,fc2)i;(-y3,fc4)> =«2^fci|fc2)(fc3|fc4). (11) 

The disorder strength is related to the mean-scattering 
time (i.e. momentum relaxation time) 

Tsm = 9 AT/c\\ ' (12) 
TTW-^A' (0) 

which in turn can be determined experimentally from 
mobility measurements using the equation 

Msm « — , (13) 

m* 

where m* is the transport effective mass and A^(0) is 
the density of states at the Fermi surface. Note that we 
have dropped a matrix-element factor in determining the 
scattering time in Eq. 1121 This is partly to account for 
multi-band effects that we are not directly accounting for 
in this paper. Since we are using mobility only as a qual- 
itative estimator of the scattering, the combination of 
Eqs. [T^ [Ql provide a good estimate of the disorder scat- 
tering w^. Although Eqs. 12 and 13 strictly apply only 
in the case of white noise disorder arising from short- 
ranged impurity scattering, we assume that they remain 
valid for long-range disorder as well in the presence of 
random charged impurity scattering since carrier screen- 
ing of the impurity potential should render the Coulomb 
disorder into an effective short-range disorder. In ad- 
dition, our use of Tsm to operationally characterize the 
semiconductor disorder should be qualitatively valid for 
any kind of impurity potential. While we have made a 
single Fermi surface approximation in Eq. SI this is more 
for the sake of simplicity. In the next paragraph, we will 
comment on how our results generalize to the multiple 
Fermi surface for example in the case of multi- 

band occupancy in the semiconductor nanowire where 
the Fermi level could lie in some high sub-band rather 
than the ground ID sub-band. 

The topological superconducting quasiparticle gap in 
the presence of a disorder potential may be determined by 
calculating the disorder averaged Green function within 
the Born approximation.^'' The disorder averaged Green 
function can be calculated from the self-energy S using 
the Dyson equation 

G~'^{k,uj)^G^\k,uj)-J:{k,uj), (14) 

where the Green function for the clean system is given 
by 

GQ^{k,uj) ^ (cj-iJo(fc)) = a;-(efe-eF,sm)T2-Afcr+-Afcr_. 

(15) 

The self-energy I](fc,aj) is approximated within the self- 
consistent Born approximation^^ as 

E„^(fc,fc') = J dki{V^x{k,ki)Vsi3{kuk))Gxs{ki)S{k-k') 

(16) 
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where {Vax{k, ki)Vsi3{ki, k)) is the disorder propagator 
from Eqs. [TUl [TTl a, /3, A, (5 are Nambu indices, and for 
brevity of notation we have dropped the w dependence 
of S and G in this equation. The self-energy within the 
first-order Born approximation is obtained by replacing 
G by Go in Eq. [TH The disorder calculation for the case 
of multiple Fermi-surfaces follows in an analogous way. 
Since momentum is a good quantum number in Eq. I14[ 
Eq. [14] does not couple the multiple Fermi-surfaces and 
therefore can be solved for each Fermi surface separately. 
The self-energy term Eq.[in]in principle has contributions 
from all Fermi surfaces. However, it is dominated by 
the contribution from the Fermi-surface with the smallest 
gap. Therefore, our calculation, which focuses on the 
Fermi surface with the smallest gap, is expected to yield 
qualitatively correct results. 

Solving the Dyson equations in the first-order Born 
approximation in the ID and 2D cases, as discussed in 
appendix \^ we obtain the disorder averaged TS quasi- 
particle gap 



(17) 



Note that the disorder renormalization of the gap van- 
ishes in the time-reversal symmetric limit since in that 
case |fc) and |— fe) are Kramer's pairs and {kp] — kp) — 0. 
The result Eq. [17] gives the reduction in the quasiparti- 
cle gap due to weak disorder within the first-order Born 
approximation. 

The Dyson equations, Eqs.[T31[Tni can in fact be solved 
for arbitrary disorder strengths within the self-consistent 
Born approximation as explained in Appendix. [B] The 
resulting TS state quasiparticle gap in the presence of 
disorder is found to be 



(18) 



where x = (ttTs^A) Kkp] ~ kp)]^^^- Note that the 
gap Eg within the first-order Born approximation in 
Eq.[T7]is not simply obtained by taking the limit of small 
X in the self-consistent Born result Eq. [THJ This is be- 
cause the Born approximation is a series expansion in 
terms of the Green function, which has a pole at uj ^ Eg, 
so that the expansion parameter in the series expansion 
diverges (at any value of x) at the relevant point and one 
cannot expect the lowest order approximation to match 
the self-consistent result. Despite, this technical point, 
the first-order and self-consistent Born approximation re- 
sults. Eg . [TTl and Eq.[TH]are qualitatively similar at small 
values of the parameter x. 

The above equation is valid as long as the expression 
under the square root is positive. In fact, the quasi- 
particle gap Eg closes for sufficiently strong disorder (or 
sufficiently small Argm) so that, by solving for Eg — 
in Eq. 1181 the disorder scattering time is found to satisfy 
the lower bound 



3.4 

Tsm > —hA^^\{kp\ - kf 



(19) 



above which the TS state quasiparticle gap Eg > 0. In 
this paper, we will also assume that Tsm > ^/ip.sm {^F,sm 
is the Fermi energy with respect to the bottom of the top- 
most filled band) need not be strictly enforced because of 
the existence of multiple bands in the actual experimental 
system, which will evade the localization problem. 

From Eqs. [13] and [H] we find that the TS state gap 
induced on the Fermi surface survives as long as the mo- 
bility of the system exceeds the threshold 



Ms 



> 3.4- 



-\{kp\ ^ kp)\^' 



(20) 



This suggests that the mobility threshold can be 
decreased by decreasing the Kramer's pair overlap 
\{kp\ — kp)\, which, as we show in Eq.[571 can be achieved 
by reducing the Zeeman potential. We will apply all these 
results to electron and hole-doped semiconductors in the 
following two sections. 



III. PROXIMITY-COUPLED ELECTRON AND 
HOLE-DOPED SEMICONDUCTORS 

In this section we review the two recently proposed sys- 
tems, electron- and hole-doped semiconductors, to which 
we will apply our results derived in the previous section. 



A. Electron-doped semiconductors 

For 2D electron-doped semiconductors with Rashba 
and Zeeman couplings, Hsm,ciean in Eq.[5]can be written 
as 



H. 



sm.cLean 



(k) 



2m* 



-a{kxcr)-z + Vza-y-ep,sm (21) 



where m* is the effective mass of the electrons, a is the 
Rashba spin-orbit coupling constant, and Vz is the Zee- 
man splitting. For a one-dimensional nanowire along the 
y-axis, the above Hamiltonian reduces to 



H sm^clean {k) 



2m* 



akax + VzOy - (.p., 



(22) 



For a chemical potential |eF,sm| < ^z, the above 
Rashba spin-orbit-coupled Hamiltonian has a single band 
at the Fermi-level with an eigenstate described by the 
spinor 



\kp) 



akp 



-^V'i+a^kl~Vz 



(23) 



Now, using Eq. [Q] the superconducting quasiparticle gap 
in the clean electron-doped case is given by 



A 



akp 



VV^+a^k'i 



:Ao. 



(24) 
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The proximity-induced pairing potential in the semicon- 
ductor, Aq, can be related to the pairing potential in the 
bulk superconductor by the relational 



Ao = A, 



A 



A + A, 



where 



\ = At\^Psc{EF) 



(25) 



(26) 



Here t represents the tunneling matrix element between 
the semiconductor and the superconductor, psc{Ep) is 
the normal state density of states of the superconduc- 
tor and A gives the rate at which electrons from the 
semiconductor tunnel into the superconductor)^ We note 
that for weak-tunneling Aq ^ A, and for strong tunnel- 
ing, Aq ~ Ag. The superconducting proximity effect in 
electron-doped InAs/superconducting Al interfaces has 
been already been observed to lead to Aq ^ A^ ^ 2 K, 
— so that one can infer that the transparency parame- 
ter A ^ As. The Kramer-pair overlap (/cf| — ^.f), which 
determines the effect of disorder scattering via Eq. [THl is 
then given by 



(27) 



which in the limit ^ 1 becomes (fc^l — kp) 



Vz 



B. Hole-doped semiconductors 

Recently, it has also been proposed^ that hole-doped 
nanowires, because of the possibly stronger spin-orbit 
coupling22ii^ and larger effective mass of holes, might 
be a better candidate for creating topological supercon- 
ductors and Majorana fermions. The hole-bands of a 
group III-V semiconductor such as InAs or GaAs are 
composed of p-orbitals, which have orbital angular mo- 
mentum L — 1. Because of strong spin-orbit coupling 
and cubic symmetry of these semiconductors, the de- 
scription of the states near the top of the valence band at 
fe = must be written in terms of spin-3/2 total angular 
momentum matrices J = L + S. The total spin-1/2 man- 
ifold of states forms a split-off band, which is separated 
by a large energy and therefore can be ignored. The re- 
sulting three-dimensional Hamiltonian is known as the 
four-band Luttinger model^, which is written as 



2rnn 



Ai + tA2)/-2A2 P»J' 



-^3 77^ PaP/sJaJ/S 

2too 



ot—x^y^z 



(28) 



where p — — iV and 71,2,3 are effective Luttinger pa- 
rameters. The effective mass mo for the holes is nega- 
tive so that the curvature of the valence bands are in- 
verted with respect to the conduction bands. We will 



consider a nanowire formed by confining the bulk hole- 
doped semiconductor confined strongly along the x and z- 
directions, so that the confinement energy scale ~ 27noL^ 
is large compared to all energy-scales we will be inter- 
ested in. Therefore one can project the three-dimensional 
Hamiltonian H30 into the lowest confinement band of the 

nanowire by replacing Px,z {px,z) = 0, p^^^ -> (pl^^) = 
2 

and adding a Rashba spin-orbit coupling proportional 
to Qi^3,3i,32 g^j^^ Zeeman term proportional to Vz to ac- 
count for electric field induced inversion symmetry break- 
ing and magnetic field-induced time-reversal symmetry 
breaking respectively. The resulting Hamiltonian for a 
hole-doped nanowire along y with a Rashba coupling (for 
the lengths in the transverse directions Lz = Lx) is given 
byr 



H 



Al 5A2 

sm,clean{k) = ^("^ + ~^ A2 J^)A;^ - akj^ 



VzJy 



(29) 



The Hamiltonian at fc = 0, com- 



where P2 

mutes with Jy. Moreover, the states with J'y — 1/4 form 
the highest energy states in the valence band. These 
states also have the lower effective mass of the two man- 
ifolds Jy = 1/4, 9/4 of states. Note that this is oppo- 
site to the two-dimensional case, where the heavy-hole 
is closer to the band-edge. Thus the light-hole states 
with Jy — ±1/2 form the set of states closest to the 
band-edge and we can restrict to this two-dimensional 
subspace Jy = ±1/2, where we identify Jy — (Jy — ±.1/2 
as a pseudo-spin variable. In this subspace, Jx ^ ctx and 
Jz ^ o'z, so that the Hamiltonian for the two upper bands 
effectively becomes 



,Ai 
2 



Hf,rr,..ri.fan{k) = —{-^ + A2)fc^ — akcF^ + VzOy + ep,, 



(30) 

The hole Hamiltonian Eq. [301 in the specific isotropic 
limit we consider, is identical (up to sign) to the 
one-dimensional electron-doped Rashba Hamiltonian in 
Eq. [21] and therefore one can expect most of the physics 
discussed in the previous sub-section for the electron- 
doped nano-wire in this limit to translate to the hole- 
doped case, except for the different effective mass and 
Rashba spin-orbit coupling. While the band-structure 
of the hole-doped wire is quite similar to the electron- 
doped case, it is not obvious from this discussion that 
the s-wave proximity effect has the same form in the 
a spin-space. In principle, the proximity-induced effec- 
tive superconducting pairing potential Aq at an inter- 
face between the semiconductor and the superconduc- 
tor can be calculated by projecting the anomalous self- 
energy Esc induced by proximity to the superconductor, 
— into the low-energy space of states \ay = ±1/2). Since 
these states form a Kramers pair, one can see that time- 
reversal symmetry only allows a conventional spin-singlet 
pairing in the ay = ±1/2 space. Furthermore, the sub- 
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space of states in the \ay = ±1/2) basis may be expanded 
in terms of the original spin eigenstates \ms = and 
\Lx,y,z = 0) (which are related to the orbital angular mo- 
mentum eigenstates jm^ — 0, ±1) by a unitary transfor- 
mation) as 



Jy - 2 



Lx — 0, Sy — - 



Ly = 0,Sy = -)] 



Jy- 2 



+ i 



Lz — 0, Sy — —- 



(31) 



Lx — 0, Sy 



(32) 



To calculate Ao, which is related to the matrix-element 
(cTy = -|-1| T,sc — —1) of a spatially smoothly vary- 
ing (on the orbital scale) proximity-induced self-energy 
one can ignore the matrix elements involving the 
orbital |La; = 0) , = 0) for SC interfaces normal to the 
l-axis. This leads to the effective barrier transparency 
parameter A in Eq. [51] to be reduced by a factor of 3. 
The corresponding pairing potential Aq from Eq. [25] is 
not significantly affected and remains the same between 
electrons and holes provided a moderately transparent 
barrier with A > As can be achieved. Therefore, more 
generally, one can expect a proximity-effect where the 
transparency-factor of the Sm/SC interface. A, is sup- 
pressed by at most a factor of order 1, which a value for 
the pairing potential Aq ~ A^, even for holes. 

To summarize, while the hole-doped nanowires have 
fundamentally different microscopic Hamiltonians con- 
sisting of a triplet of orbitals with p-symmetry, for 
nanowires it is possible to work in the limit of symmet- 
ric confinement along the two directions transverse to 
the wire, where we have shown that the effective Hamil- 
tonian and proximity-effect are identical in form to the 
electron-doped wire case. Of course, the parameters in 
the Hamiltonian are expected to be quite different and in 
particular the effective mass and Rashba spin-orbit cou- 
pling are expected to be much larger. At the same 
time, the corresponding mobility for the semiconductors 
that can be achieved are found to be significantly lower, 
— indicating much stronger disorder for hole-doped sys- 
tems. Thus, the hole systems have the advantage (disad- 
vantage) of stronger SO coupling (larger disorder), and 
whether the net prospect for realizing the non-Abelian 
Majorana mode is greater in hole-based systems (than 
the electron-based systems) or not will depend entirely 
on which of these two effects wins over in realistic sam- 
ples. Another possible advantage of the hole systems is 
the higher value of the Lande g-factor which necessitates 
weaker magnetic field values in achieving the condition 



for topological superconductivity. One should, however, 
realize that very little actual experimental work exists in 
hole-based semiconductor nanowires. However, it is en- 
couraging that the superconducting proximity effects in 
both electron and hole-based semiconductors have been 
observed in experiments i^i^ Therefore, we do believe, 
that the hole-based systems (e.g. InAs, InSb, SiGe, and 
perhaps even GaAs) should be included in the experi- 
mental search for the Majorana mode in semiconductor- 
superconductor heterostructures . 



IV. TS STATE QUASIPARTICLE GAP IN THE 
PRESENCE OF DISORDER 



The clean-limit quasiparticle gap A in the topological 
superconducting state can be obtained for the Rashba 
spin-orbit coupled Hamiltonian Eq.[5T] using Egs.lMland 
[21] as 



A = 



AA, 



aki. 



A+A, yrf^ 



(33) 



This expression suggests that to obtain a large gap, A, 
in the semiconductor one should restrict oneself to the 
limit where akp ^ Vz- Moreover, to remain in the req- 
uisite odd sub-band (confinement bands) limit it is also 
necessary to choose a chemical potential ep^sm in the gap 
(at fc = 0) induced by the Zeeman splitting Vz^ This 



restricts kp 



in the highest filled sub-band. In this 



limit one can approximate akp ~ Eso = 2m* a^, which 
we call the spin-orbit energy Eso ■ This spin-orbit energy 
Eso is a crucial energy scale in the problem which is to 
be compared with the other independent energy scales 
such as Vz, As and l/rsm- Therefore, from Eq. 30, to 
obtain a large quasiparticle gap at the Fermi surface the 
system must be restricted to the regime 



yz J5 Esq. 



(34) 



The semiconductor system is in the topological state if 
and only if the Pfaffian invariant of the relevant BdG 
Hamiltonian is negative)^ which is equivalent to the 
topological condition 



Vi > A^ + g2 , 



(35) 



where Vz = Vz/(1 + A/A^), Ao = Ao and ep^sm = 
/(I + A/ As) are the renormalized parameters in 
the effective Hamiltonian of the proximity coupled semi- 
conductor after the superconductor has been integrated 
out4i21 Thus, using Eq. [25] and Eq. [35] the topological 
condition in terms of the bare parameters can be written 
as. 



(36) 



Eq. [3i and Eq. [Ml together with Eq. [18] (or Eq. [20]), 
define the optimal constraints on the material considera- 
tions for the realization of the TS phase in the semicon- 
ductor. 
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In what follows, we assume that the system is rea- 
sonably deep in the topological superconducting phase 
(i.e. Vz > 2A, 2ei?,sm) so the smallest gap occurs near 
the Fermi wave- vector |A;| ^ kp instead of near |A;| ~ Oi^ 
Therefore in attempting to optimize the quasiparticle gap 
one notices from Eqs. [SH [36]a hierarchy of energy scales 



Eso >Vz >2X>2A> 2Eg 



(37) 



The first two-inequalities from the left are constraints 
between the external parameters i?so: and A, out of 
which Vz and A must be tuned to satisfy the above in- 
equality. The other external parameter is the disorder 
scattering rate in the semiconductor, t~^, which must 



1 



The two- 



be reduced to satisfy A > (nTsm) 

inequalities on the right-hand side of Eq. 1371 are between 
the derived quantities A and Eg, which are completely 
determined once Eso,Vz , X,Tsm are given, so that the 
condition 2A > 2A and 2A > 2Eg follow as direct conse- 
quences of Eqs. 133] and Eq. [151 respectively. Therefore, 
the disordered TS state quasiparticle gap, Eg, is con- 
strained by the spin-orbit energy scale EgQ, which must 
be as large as possible in the problem. 

In the presence of disorder, one can use Ea.lT^together 
with Eq. [27] to show that for Vz ^ Eso , the disorder 
renormalized TS state quasiparticle gap is given by 



E„ = AVi -3x + 3x2 _ 32,3 



(38) 



-2/3 



4/3 



The mobility 



where x = (ttt^^A) . , 

V V ^z+-'^so , 

threshold, which is the minimum mobility required to 
realize a non-zero TS state gap Eg in the presence of 
disorder, is obtained by using Eqs. [501 [IZ] to be 



fism > 6.8- 



ehVzjAs + Vz/2) 
^*AsEsoVvf 



^SO 



(39) 



which vanishes in the limit Vz 0- Note we have rein- 
troduced the ?i = 1 to make the expression dimensionally 
consistent. Since fism increases with A, which in turn 
increases with A, for deriving Eq. 1391 we have assumed 
the A is chosen to reach its maximum value A = Vz/2, 
above which the system can no longer be considered to 
be deep in the topological phase (see Eq.[271). In fact, we 
will assume that Vz has been tuned to be Vz = 2\ for 
the rest of the paper. 

Of course, from Eq.[3ni there is a threshold value for Vz 
only above which the TS state itself is realized. There- 
fore, Vz cannot be taken as zero in Eq. [35] to maximize 
the disordered TS state gap Eg. In Eq. [38] to deduce 
the optimum Eg, we use the r.h.s of Eq. [33] to substi- 
tute for A and, for a robust topological state, we take 
A = Vz/2. The r.h.s of Eq. [25]is then a function of Vz, 
Tsm (or fism ~ GTsm/m*), Aj and Eso- Among these 
A^ (the gap in the adjacent superconductor) and Eso 
(the spin-orbit coupling energy in the semiconductor) are 
known experimental quantities. Therefore Eq. 1381 defines 
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FIG. 1: (a) Calculated disordered TS state quasiparticle gap 
Eg as a function of semiconductor mobility /ism for electron- 
doped wires with Eso = 2m* — 2K^ and Eso = 4:K al- 
lowing for future electric-field enhanced spin-orbit coupling.— 
(b) The values of Vz = 2A realizing Eg in panel (a) plotted 
with Ham- The gaps in the adjacent bulk superconductor are 
taken to be Aa = 2K corresponding to Al and As — 
corresponding to Nb. 



a formula for Eg in terms of Vz parametrically related 
by fism ■ For given values of figm , we maximize Eg with 
respect to Vz and plot these maximized values of Eg as a 
function of fism for electron and hole-doped semiconduc- 
tors in Figs. [1] and [21 respectively. These plots arc the 
central results of this paper. 

From Fig. [T] it is clear that for electron-doped wires 
with a mobility of 100,000 cm^/V-s^^, which has been 
achieved for wires with a diameter 60 nm, a measurable 
robust TS state gap of 50 — 100 mK is achievable. As is 
also clear from Fig. [21 for hole-doped wires with a signif- 
icantly larger spin-orbit energy Eso ^ 30 meV, a larger 
gap of 0.8 K is achievable for a mobility of 100 cm^/V- 
It should be noted that in general the holes have 
much lower mobilities than electrons in the same semi- 
conductor material even after taking into account their 
large effective mass difference. 



V. SUMMARY AND CONCLUSION 

In this paper we have addressed the effects of dis- 
order from the semiconductor surface on the stability 
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FIG. 2: (a) Calculated disordered TS state quasiparticle gap 
Eg as a function of semiconductor mobility for hole-doped 
wires with Eso = 300KM->^ (b) The values of Vz = 2A re- 
alizing Eg in panel (a) plotted with /ism- The gaps in the 
adjacent bulk superconductor are taken to be As — 2K cor- 
responding to Al and As = iK corresponding to Nb. 



of the topological superconducting state on ID and 2D 
electron- and hole-doped semiconductors proximity cou- 
pled to bulk 3D s-wave superconductors. In recent 
works^'' it has been shown that the effects of disorder 
from the adjacent bulk superconductor on the TS state 
is minimal, so we have ignored this effect. However, since 
the TS state in the semiconductor explicitly breaks the 
time-reversal symmetry (due to an external Zeeman split- 
ting), even the disorder from the semiconductor surface 
itself has an effect on the stability of the TS state and will 
close the BCS-like superconducting quasiparticle gap at 
the semiconductor Fermi surface for low enough mobility. 
We have shown that the disordered TS state quasiparti- 
cle gap Eg is suppressed with increasing values of the 
ratio r between the Zeeman splitting Vz and the spin- 
orbit energy scale Eso = 2m* a^, r = Vz/Eso- The 
dimensionless quantity r measures the extent of the TR 
symmetry breaking in the semiconductor to produce the 
TS state and Majorana fermions. 

Our main results are plotted in Figs. [T] and [2] They 
display the achievable disorder renormalized TS state gap 
Eg as a function of the semiconductor mobility fj,sm in 
both electron and hole-doped semiconductor wires. From 
Fig. [1] we find that for electron-doped wires with a mo- 
bility of 100, 000 cm^/V-s"^^, which has been achieved for 



wires with a diameter 60 nm, a measurable robust TS 
state gap of 50 — 100 niK is achievable. From Fig. [5] we 
find that for hole-doped wires with a significantly larger 
spin-orbit energy Eso ^ 30 meV, a larger gap of 0.8 
K is achievable for a mobility of 100 cm^/V-Si^ Thus, 
despite the fact that the semiconductor TS state breaks 
the time-reversal symmetry, realizing a robust TS state 
is possible in both electron- and hole-doped systems even 
in the presence of realistic disorder on the surface of the 
semiconductor. 

Before concluding, it may be useful to provide a qual- 
itative discussion and a general perspective on the is- 
sue of optimal materials for the realization of the pre- 
dicted TS phase and Majorana modes in semiconductor- 
superconductor sandwich structures. The first point to 
emphasize is that there are far too many independent 
physical parameters entering the problem for theory to 
be particularly useful in this context. At least six inde- 
pendent material parameters play a role here: a, kp, Vz, 
eF,sm, As, t (which determines A in Eq.[55]). In addition, 
disorder, which plays a crucial role by the non-existence 
of Anderson's theorem in the time-reversal invariance 
broken (i.e. Vz ^ 0) situation, is itself characterized 
by several independent physical mechanisms: disorder in 
the adjacent bulk superconductor, long-ranged disorder 
in the semiconductor, disorder-induced fluctuations in t, 
interface disorder, puddles, and so on. A physical system 
characterized by such a large (> 10) set of independent 
parameters, all of which are important in determining the 
fate of the TS state, is not an easy optimization problem 
to deal with theoretically. In particular, it may turn out 
to be more practical to simply carry out experiments on 
available systems (e.g. InAs or InSb on Al or Nb) than 
to trust precise theoretical quantitative predictions with 
respect to the optimal materials choice. 

With these important caveats in mind, we can, 
nonetheless, make some general qualitative statements, 
quite apart from the more precise quantitative state- 
ments made in the main body of this work. In the clean 
case, where disorder is not an issue (i.e. extreme high- 
mobility semiconductor systems), it is obvious that the 
semiconductor spin-orbit coupling should be as large as 
possible. This continues to be the case even in the pres- 
ence of disorder, making Esq (« akp) to be the key phys- 
ical parameter, which should be made as large as possible 
(including enhancement by external gates if feasible — ) 
with no constraints. Similarly, disorder is always harm- 
ful to the TS state, and therefore 1 /vsm should be made 
as small as possible by working with the highest mo- 
bility materials. Although having a spin-orbit coupling 
strength as large as possible is obviously helpful to the re- 
alization of the TS phase, this may have the detrimental 
effect of enhancing impurity-induced scattering between 
the spin-split bands, thus reducing the mobility)^ and 
therefore some optimization may be necessary even with 
respect to spin-orbit coupling and mobility in order to 
satisfy the general requirement of having both the high- 
est possible spin-orbit coupling and mobility. 
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The same is to a limited degree true of the supercon- 
ducting gap As (or equivalently T^) in the bulk superc- 
noductor, but this advantage is considerably nullified by 
the fact that the proximity effect is determined primar- 
ily by the parameter A (see Eq. 1251) , and this parameter 
should be made as large as possible, even at the cost of 
having a lower in the bulk superconductor (e.g. Al 
may be preferable over Nb if Aai > Awb)- We note that 
having an arbitrarily large proximity gap Aq in the semi- 
conductor may not necessarily help unless Eso can also 
be simultaneously increased (e.g. by an external gate) 
since Eso > Vz > 2X > 2A > 2Eg is a necessary condi- 
tion for the existence of the TS phase. Similarly, having 
a very small Vz does not help since Vz > 2A determines 
the condition for a robust topological phase, and having 
a small Vz necessarily requires having a much smaller A, 
which severely restricts the temperature range where the 
experiments investigating the topological phase can oc- 
cur since T < A is necessary. For the same reason a small 
Eso, e.g. GaAs, with correspondingly smaller Vz and A 
would not work either since the operational temperature 
would then be too small. 

As for disorder, our discussion in this work has con- 
sidered the semiconductor disorder to be short-ranged, 
which is consistent with interface disorder (or roughness) 
scattering which is known to be important in semicon- 
ductor nanowires. Carrier screening should render long- 
range Coulomb disorder arising from random charged im- 
purities in the semiconductor into effective short-ranged 
disorder, making our model of rather general validity. 
The fact that we model disorder through an effective re- 
laxation time Tsm defined through the measured carrier 
mobility makes the model robust with respect to differ- 
ent types of disorder in the semiconductor, but we note 
that for long-ranged Coulomb disorder is really only 
a lower bound on the effective disorder. (Disorder arising 
from impurities in the superconductor does not affect the 
TS phase at all, see Ref.~ for details.) 

A simple dimensional argument, which has been au- 
thenticated by a renormalization group calculation;^ 
suggests Atsju > 1, i-e. t~J^ < A, as the bound on 
the effective scattering rate for the topological gap A to 
survive disorder in a spinless chiral p-wave superconduc- 
tor characterized only by these two parameters (unlike 
our problem which is characterized by many physical pa- 
rameters). We note that our microscopic Born approx- 
imation calculation (Eq. [T9| is consistent with this case 
deep in the TS phase for Vz ^ Eso- For Eso > Vz, 
the effect of disorder scattering is suppressed so that the 
topological gap survives disorder for a broader range of 
parameters Ar^m > {Vz / Eso)'^ ■ Unfortunately for the 
semiconductor-superconductor sandwich structures of in- 
terest to us, ATsm > (Vz/Eso)'^ only provides a neces- 
sary condition since A cannot be made arbitrarily large 
unless Eso i-e. the spin-orbit coupling is itself arbitrarily 
large. The sufficient conditions for the gapped topolog- 
ical state to exist in the semiconductor sandwich struc- 
tures are given by Eqs. 33, 34, which require "large" Vz 



making it necessary to have even larger Eso^- Thus, the 
existence of the TS phase in the realistic semiconductor 
sandwich structures depends on the seamless and multi- 
parameter materials optimization of spin-orbit coupling, 
spin splitting, superconducting gap, and disorder in an 
inter-dependent manner. 

We conclude with the observation that the single 
most important ingredient for the TS state to emerge 
in semiconductor-superconductor sandwich structures is 
to have a large spin-orbit coupling and high mobility. 
For electrons both InAs and InSb satisfy the necessary 
requirements if the mobility can be enhanced around 10^ 
cm^/V-s. For holes the spin orbit coupling being very 
large, much lower mobilities (^ 10^ cm^/V-s) should 
work. These considerations imply that holes should be 
looked into experimentally for the realization of the TS 
phase. In addition, having several (odd number of) 
subbands^ occupied in the nanowire should help since 
it reduces disorder through screening and enhances the 
effective values of Eso- We emphasize that the exis- 
tence of the TS state depends delicately on the competi- 
tion among superconducting proximity effect, spin-orbit 
coupling, spin splitting, and chemical potential. (Disor- 
der scattering always suppresses the TS state.) As such, 
any materials optimization must carefully choose among 
these physical parameters while at the same time main- 
taining the highest possible mobility (i.e. smallest possi- 
ble disorder scattering) . The experimental search for the 
Majorana mode should thus concentrate on semiconduc- 
tor systems of the highest mobility and the strongest SO 
coupling using the appropriate superconductors which 
produce the best quality interface between the supercon- 
ductor and the semiconductor. 
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Appendix A: Dyson's equation and disorder 
averaged Green's function 

In this appendix, we start by first simplifying the 
Dyson equations Eqs. [14] and [16] and then solve the 
equations for the first order Born approximation case. 
The self-consistent Born approximation is solved in Ap- 
pendix. [B] Substituting Eqs. [TU] and [11] into Eq. [TH] the 
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self-energy E is written as 



where 



S(fc) 



■G(fci) 



2 



C= de 



(A8) 



(fc|fci)*(-fc|-fei)*il-^G(fci; 



2 

(1 + r, 



(Al) 



and az{k,uj) — (e^ ) remains unrenormalized. 

The lowest-frequency pole in the disorder averaged 
Green function G{k,u!) given by Eq. IA41 which corre- 
sponds to the quasiparticle gap, occurs when ek — ep.sm — 

= and ag — \a-^-\■ The disorder averaged quasiparti- 
cle gap io can then be obtained by solving 



In principle, the disorder averaged scattering po- 
tential {Va\{k,ki)Vsp{ki,k)) involves integrals of 4 
Bloch eigenstates {r,a\k) — Uk{r;a) of the form 
v'^ J drJ2^yul,{r;a)uk{r;a)ul{r;a')uk'{r;a') where 
the sum over a implies that the disorder is spatially 
local spin-independent disorder. Having several sub- 
bands occupied in the ID nanowire, as has recently 
been propose d^ ''1'^° , should help the situation consider- 
ably, both by expanding the regime of stability of the TS 
phase and by screening out the disorder due to higher 
carrier density. Substituting the wave-function into the 
matrix elements, we get 



LU[1 



fo 



VA2 - a;2 



Afl 



fi 



(A9) 



For sufficiently weak disorder r^^^ ^ A, the quasiparticle 
gap cj ^ A, so that we can expand in the reduction of the 
quasiparticle gap x = A — a; to lowest order and obtain 
the quasiparticle gap reduction to be 

-/ ^-l\2/3 



A 

21/3 



(,/o - /l) 



2/3 



(AlO) 



{\k\,9k\\k'\,ek')'=^{kF,9k\kF,ek' 
= (0| e*(^--^)(^'=-^'=') |0) = fiOk - 



This leads to the pole of the disordered average Green 
function G which is also the TS state quasiparticle gap 
in Eq. [JT] 



(A2) 



1 



l^-lj^inT,„,A)-^/^\{kF\~kF)\ 



|4/3 



(All) 



Substituting the matrix elements into the disorder- 
induced self-energy we find that. 



E(fc,w) 



dk,[\f{ek~0k^)\ 

(1 + ^^)^.. ^(l 



2(G(fcl)-f T,G(fcl)T,) 



fi9k-ek,f^^^G{k,) 



fi&k — Oki) 



-G(fci)- 



(A3) 



Appendix B: Analytic solution of the self-consistent 
Born approximation 

In this appendix we provide an exact solution for the 
quasiparticle gap within the self-consistent Born approx- 
imation. To do this, we start with the self-consistent 
version of Eq. IA51 



ao{k,uj) 



To solve Dyson's equations, it is convenient to define the 
functions aQ^z,+.-ik,uj) by the relation 

G^^{k,uj) — aQ{k,uj)+az{k,LL!)T2+a^{k,uj)T^+a_{k,uj)T_. — foiT^Tsm)^^ / de' 

(A4) 

Substituting this into the Dyson equations (Eq. [T4|) to- 
gether with the expression for the self-energy (Eq. IA3|) in 
the lowest order Born approximation we get 



ao(fc, Lo) 



ao(fe,iLj)2 - (e' - a^(A;,w))2 - \a+{k,uj)\^^ 
az{k,uj) = (e - ep^sm) 

(e' — CF.sm) 



ao(fc,a;) = a; - /o(7rTs„) ^ / de' 



a+ik,Lu) = Ae'"'' 



de' 



ao(fc,w)2 - (e' - a4fe,a;))2 - |a+(fc,w)|2 



u;2 _ (e' _ ep^s^r - |A|2 



az{k,uj) = {e - ep^sm) + fo{T^Tsm) ^ de' 



a+{k,uj) 



{^'-^F,sm) V^''ao(fe,c^)2-(e'-a.(fe,c^))2-|a+(fe,c^)|2 



oj^ - {e' - ep^sm)^ - \A\^ 



(Bl) 



a+(fc,w) = Ae^"" - fliTTTsmY 



de'- 



As '° Performing the e' integral in the second equation we im- 

(e' — ep^smY — |A|fiediately obtain az{k,uj) = 0. Performing the e' integral 



(A5) 



where /o = / f |/(0)P, /i = / §fi9)^e'''. The integrals 
in Eq. IA5I can be written analytically as 

ao = w + fainTsmy^C^ (A6) 
a+ = [A + h{nTs^)-\A]e'''' (A7) 



in the other 2 equations 

ao — iu + U' 



ao 



a+ 



A- 



(B2) 
(B3) 
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where the square root is chosen so that the real part 
of \Ja\- al, a = TrfoinTsm)'^ and b = 7r/i(7rrs„)"^ 

Writing y — ^ a\ — a^, the pair of equations can be com- 
bined to a single equation for y 



1 = 



A2 



{y - hY {y - af 



(B4) 



At cj ^ 0, the equation has 4 real roots y a and y ~ 
6 ± A. The quasiparticle gap is defined by the value 
of oj above which the Green function G{uj) has a finite 
imaginary part. This corresponds to the value of w where 
2 of the real roots of y merge into a pair of complex 
roots, so that Eq. IB4I has 2 real and 2 complex roots. 
Using the properties of quartic polynomial equations, the 



transition point between 4 real roots and 2 real roots, 
is signalled by the vanishing of the discriminant of the 
quartic polynomial. Solving for the frequency lo where 
the discriminant vanishes leads to the solution 

(B5) 

where g — — . In the 2D Rashba and ID supercon- 
ducting case, we get the solution of the TS state quasi- 
particle gap in Eq. [18] 

= A[l - 3 (^T,™A)-'/^ - kF)f^ 

+ 3 (ttt.^A)-*/^ \{kF\~ kF)\^'^ - 3 (ttt.^A)-" \{kp\- kF)\^]"^. 

(B6) 
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